ESC-TR-99-066 


Technical  Report 
1063 


Calculation  of  Resident  Space  Object  Color 
Temperature  and  Emis  s ivity- Ar  ea  from 
MSX  SPIRIT  in  Infrared  Data: 
Emissive  Reference  Sphere  Results 


R.L.  Lambour 


13  September  2002 


Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
Lexusgtos ,  Massachusetts 


Prepared  for  the  Department  of  the  Air  Force  under  Contract  FI 9628-00- C-0002. 
Approved  for  public  release;  distribution  is  unlimited. 


A- DA 


This  report  is  based  on  studies  performed  at  Lincoln  Laboratory,  a  center 
for  research  operated  by  Massachusetts  Institute  of  Technology  .  This  work 
was  sponsored  by  the  Department  of  the  Air  Force,  SMC,  under  Contract 
F19628-00-C-0002. 

This  report  may  be  reproduced  to  satisfy  needs  of  U.S.  Government  agencies. 


The  ESC  Public  Affairs  Office  has  reviewed  this  report,  and 
it  is  releasable  to  the  National  Technical  Information  Service, 
where  it  will  be  available  to  the  general  public,  including 
foreign  nationals. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 

FOR  THE  COMMANDER 


Gan^  I'utungian 
Administrative  Contracting  Officer 
Plafts  and  Programs  Directorate 
Contracted  Support  Management 


Non~lincoln  Recipients 

PLEASE  DO  NOT  RETURN 

Permission  is  given  to  destroy  this  document 
when  it  is  no  longer  needed. 


Massachusetts  Institute  of  Technology 
Lincoln  Laboratory 


Calculation  of  Resident  Space  Object 
Color  Temperature  and  Emissivity-Area  from 
MSX  SPIRIT  III  Infrared  Data: 
Emissive  Reference  Sphere  Results 


R.L.  Lambour 
Group  91 


Technical  Report  1063 


13  September  2002 


Approved  for  public  relase;  distribution  is  unlimited. 


Lexington 


Massachusetts 


ABSTRACT 


This  report  discusses  the  calculation  of  color  temperature  and  emissivity-area  product  from 
resident  space  object  (RSO)  observations  that  were  obtained  during  a  series  of  surveillance 
experiments  carried  out  with  the  Midcourse  Space  Experiment  (MSX)  spacecraft.  During  these 
experiments,  simultaneous  infrared  and  visible-band  data  were  collected  with  the  SPIRIT  III 
LWIR  infrared  radiometer  and  the  Space-Based  Visible  (SBV)  instruments,  respectively. 
Algorithms  were  developed  for  the  SPIRIT  III  Data  Reduction  pipeline  that  allow  calculation  of 
color  temperature  and  emissivity-area  product  for  RSOs  given  that  the  RSO  was  observed  in  at 
least  two  wavelength  bands.  This  report  discusses  the  algorithms  and  characterizes  their  accuracy 
using  data  taken  on  five  emissive  reference  spheres  that  were  deployed  from  the  satellite  over  the 
lifetime  of  the  SPIRIT  III  instrument. 
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1.  INTRODUCTION 


The  Midcourse  Space  Experiment  (MSX)  spacecraft  was  launched  in  April  1996  into  a  circular, 
898-km,  sun-synchronous  orbit  with  an  inclination  of  99  deg.  The  primary  purpose  of  the  satellite  was  to 
collect  target  and  background  phenomenology  data  in  order  to  facilitate  detection,  acquisition,  and 
tracking  of  ballistic  missile  targets  in  midcourse  flight,  as  well  as  discrimination  of  lethal  and  non-lethal 
ballistic  missile  targets.  The  spacecraft  is  equipped  with  a  suite  of  optical  instruments  that  make 
measurements  from  the  far  ultraviolet  (1 10  nm)  to  the  long  wave  infrared  (26  pm)  region  of  the  spectrum. 
The  optical  instruments  are  the  SPatial  InfraRed  and  Imaging  Telescope  (SPIRIT)  III,  a  cryogenically- 
cooled  radiometer  and  interferometer,  the  Space-Based  Visible  (SBV)  camera,  and  the  Ultraviolet/Visible 
Imaging  and  Spectrographic  Imaging  (UVISI)  sensors. 

During  the  10-month  period  before  the  SPIRIT  III  cryogen  was  exhausted  (April  1996  through 
February  1997),  use  of  the  spacecraft  was  shared  between  eight  Principal  Investigator  (PI)  groups.  This 
time  period  is  referred  to  as  MSX  Phase  I  data  collection  or  simply  as  Phase  I.  The  Surveillance  PI  group 
executed  fifteen  experiments  (data  collection  events  or  DCEs)  during  Phase  I.  The  experiments  utilized 
SPIRIT  III  and  SBV  simultaneously  and  are  summarized  in  Reference  1.  Five  different  types  of 
experiments  were  run,  all  of  which  collected  data  simultaneously  in  the  infrared  (IR)  and  visible  bands. 

The  primary  objective  of  the  experiments  was  to  examine  the  multi-band  IR  and  visible  data  for 
its  utility  in  supporting  Space  Object  Identification  (SOI)  and  space  surveillance.  The  objective  of  SOI  is 
draw  conclusions  about  RSOs  (what  it  is  and  what  it  is  doing)  using  the  observed  and  inferred 
characteristics.  Those  characteristics  include  the  RSOs  temperature,  area,  emissivity,  absorptivity, 
reflectivity,  and  temporal  trends.  The  Lincoln  Laboratory  SPIRIT  III  Data  Reduction  pipeline1  calculates 
two-color  target  temperature  and  emissivity-area  product  (eA)  routinely  for  all  RSO  observations  that 
have  been  correlated  to  known  objects.  This  report  describes  the  algorithms  utilized  by  the  pipeline  to 
calculate  T  and  eA  and  discusses  the  accuracy  to  which  they  can  be  detrmined.  Section  2  describes  the 
SPIRIT  III  radiometer.  Section  3  presents  the  physics  behind  the  determination  of  an  RSOs  temperature 
and  its  IR  signature.  Section  4  presents  the  two-color  temperature  calculation  algorithm  and  Section  5 
discusses  the  accuracy  to  which  it  can  be  determined.  Section  6  describes  the  calculation  and  accuracy  of 
eA.  Section  7  presents  a  discussion  of  the  results  and  Section  8  presents  the  conclusions  of  this  report. 
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2.  THE  SPIRIT  III  RADIOMETER 


The  SPIRIT  III  sensor  is  the  primary  sensor  on  MSX.  It  was  designed  and  built  at  the  Utah  State 
University  Space  Dynamics  Laboratory  in  Logan,  Utah.  The  telescope  covers  the  spectrum  from  the  mid¬ 
wave  IR  (MWIR)  to  the  long-wave  infrared  (LWIR).  SPIRIT  III  consists  of  an  off-axis,  high  stray-light 
rejection,  reimaging  telescope  with  a  37-cm  diameter  aperture,  a  6-band  scanning  radiometer,  and  a  6- 
channel,  high-resolution  interferometer-spectrometer,  all  of  which  are  cooled  to  cryogenic  temperatures 
(~1 0-13  °K)  by  a  dewar/heat  exchanger  which  was  filled  with  solid  hydrogen.  The  cryogen  lasted  10 
months  once  the  spacecraft  was  on  orbit. 

The  interferometer  was  not  used  during  the  space  surveillance  experiments  and  will  not  be 
discussed.  The  SPIRIT  III  radiometer  is  a  line-scanning  radiometer  that  measures  IR  radiation  in  6 
spectral  bands  that  range  from  4.2-26  |im  and  have  a  spatial  resolution  of  90  |irad.  The  wavelength 
bands  are  listed  in  Table  1.  The  wavelength  bands  that  were  processed  for  the  surveillance  experiments 
are  shown  in  bold  type.  A  bi-directional  scan  mirror  is  used  to  build  up  an  image  over  a  variable  field  of 
regard  that  is  centered  on  the  telescope  boresight.  The  radiometer  scan  mirror  can  remain  fixed  or  can 
operate  at  a  constant  scan  rate  of  0.46  deg/sec  with  a  selectable  field  of  regard  of  1  x  0.75°,  1  x  1.5°,  or  1 
x  3°.  Therefore,  the  time  required  to  build  up  a  scan  varies  from  1.63  sec  to  6.52  sec,  depending  on  the 
scan  mode  selected.  The  scan  mode  used  during  the  surveillance  experiments  generated  1  x  1 .5°  scans  at 
a  “frame  rate”  of  1  scan  every  3.26  sec.  These  individual  scans  are  referred  to  as  scenes. 

The  integration  time  (the  time  required  to  read  out  all  192  rows  of  the  array  columns)  was  also 
selectable.  The  surveillance  experiments  used  a  mirror  scan  mode  with  an  integration  time  of  ~2.7  msec. 
Given  the  scan  rate  of  the  mirror,  this  produced  an  effective  dwell  time  on  a  pixel  of  ~1 1  msec.  Thus,  this 
mode  of  operation  over-sampled  the  field  of  regard  by  about  a  factor  of  4  (~4  samples  were  collected  for 
any  point  in  the  field  of  regard).  This  allowed  the  instrument  to  build  up  enough  SNR  to  detect  relatively 
faint  objects.  The  sensitivity  of  each  radiometer  band  is  characterized  in  Table  1  by  the  noise  equivalent 
irradiance  (NEI).  Additional  details  on  the  radiometer  can  be  found  in  References  2  and  3. 

The  radiometer  underwent  extensive  ground  and  on-orbit  calibration  in  order  to  determine 
calibration  quantities  (e.g.,  NEI,  responsivity,  radiometer  pixel  linearity,  etc.).  It  is  beyond  the  scope  of 
this  report  to  describe  the  calibration  procedures  and  activities.  For  additional  information,  the  reader  is 
referred  to  References  4  and  5  and  the  references  therein. 

TABLE  1 

SPIRIT  III  Band  Characteristics  (FWTM§) 


Spectral  Band 

Wavelength  (pm) 

NEI  (W/cm2)* 

A 

6.03-10.92 

3.2x1  O'18 

B1 

4.21-4.37 

4.5  x  10'17 

B2 

4.23-4.47 

3.7  x  10‘17 

C 

11.05-13.27 

1.4  x  10‘17 

D 

13.32-16.05 

5.4x1  O'18 

E 

18.03-26.79 

1 .5  x  10‘17 

*FP  Temperature  =  12  K 
§Full-width  tenth  maximum 
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3.  RSO  TEMPERATURE  AND  IR  SIGNATURE 


The  temperature  of  an  RSO  is  influenced  by  a  number  of  factors  that  are  diagrammed  in  Figure  1 . 
Consider  the  concept  of  instantaneous  thermal  equilibrium,  in  which  the  amount  of  energy  emitted  by  the 
RSO  is,  at  any  instant,  equal  to  the  amount  of  energy  absorbed.  This  situation  represents  the  conservation 
of  energy.  The  power  (energy  per  unit  time)  incident  on  the  surface  of  the  RSO  can  be  written  as: 


Incident  Power  =  Ap&sun  +  aFApQ>sun  +  FAp®E*nh  +  Qu 


(1) 


where  F  is  the  Earth  view  factor,  Ap  is  the  projected  area  of  the  RSO  (e.g.,  kx2  for  a  sphere),  Osun  is  the 
energy  flux  incident  on  the  RSO  from  the  Sun  (-1360  W/m2  at  the  Earth),  d>Earth  is  the  energy  flux 
incident  on  the  RSO  due  to  the  warm  Earth,  a  is  the  albedo  of  the  Earth  scene  that  is  visible  to  the  RSO, 
and  Qmt  is  a  term  describing  the  internal  heat  generation  of  the  RSO.  The  Earth  view  factor  is  the  fraction 
of  4k  solid  angle  subtended  by  the  Earth  at  the  location  of  the  RSO.  The  first  term  on  the  right  of  the 
equal  sign  represents  the  power  incident  from  direct  solar  illumination.  The  second  term  represents  the 
power  incident  on  the  RSO  due  to  reflection  of  solar  illumination  by  the  Earth.  The  albedo,  ay  represents 
the  reflectivity  of  the  Earth  averaged  over  the  entire  Earth-scene  visible  to  the  RSO.  The  albedo  is  highly 
variable,  depending  on  the  visible  terrain  type  (e.g.,  snow,  desert,  forest,  ocean,  etc.)  and  the  amount  of 
cloud  cover.  The  Earth’s  albedo  can  range  from  0.08  to  0.8,  although  0.3  is  an  accepted  average  value6. 
The  third  term  represents  the  contribution  of  upwelling,  terrestrial  thermal  (long-wave  IR)  radiation 
(Earthshine)  incident  on  the  RSO  due  to  its  proximity  to  the  warm  Earth.  This  upwelling  power  is 
principally  dependent  on  the  temperature  of  the  Earth’s  surface  and  atmosphere  and  the  amount  of  cloud 
cover  .  The  final  term  represents  the  contribution  of  internal  heat  sources.  The  power  absorbed  by  the 
RSO  surface  is  a  fraction  of  the  incident  power: 


Absorbed  Power  =  A p$>sw  +  aFApOsun  +  gFAp&Eanh)  +  CkQmt  (2) 

where  Os  represents  the  solar  absorptivity  coefficient  for  the  RSO  surface  material.  In  general,  this 
coefficient  varies  with  wavelength;  the  relation  in  Equation  (2)  assumes  that  the  RSO  is  a  graybody.  A 
graybody  is  an  imperfect  Planckian  radiator  for  which  the  absorptivity  and  emissivity  are  independent  of 
wavelength.  The  graybody  absorptivity  and  emissivity  values  are  also  less  than  1 .  We  have  assumed  in 
Equation  (2)  that  the  absorptivity  of  the  interior  surface,  (%,  may  be  different  from  the  solar  absorptivity  of 
the  exterior  surface.  Equation  two  also  includes  a  factor,  g,  where  go^  =  (%„  the  absorptivity  at  infrared 
wavelengths  (where  0Earth  is  largest).  For  a  graybody,  we  assume  g  -  1 . 
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Earth 


Figure  1.  Diagram  of  RSO  thermal  environment.  The  energy  flux  from  each  source  is  shown,  as  well  as  the  energy 

flux  from  the  radiating  RSO. 


The  power  emitted  from  the  RSO  surface  can  be  approximated  as: 

Emitted  Power  =  A 5d>Rso  -  QrAsCfT1  (3) 

where  As  is  the  total  surface  area  of  the  RSO  (e.g.,  4jtr2  for  a  sphere),  €,r  is  the  infrared  emissivity  of  the 
RSO  surface  material,  <7 is  the  Stefan-Boltzmann  constant  (5.67  x  10'8  J/m2-K4),  and  T  is  the  temperature 
of  the  RSO.  In  general,  the  emissivity  varies  with  wavelength,  but  Equation  (3)  assumes  that  the  RSO  is 
a  graybody. 

Assuming  thermal  equilibrium,  Equations  (2)  and  (3)  can  be  equated.  Solving  for  the  temperature 
shows  that  the  equilibrium  temperature  of  an  RSO  depends  upon  its  cross-sectional  area,  and  the  optical 
properties  of  its  surface  material. 

The  temperature  of  the  RSO  is  estimated  by  examination  of  the  RSO  signature  observed  by  the 
SPIRIT  III  radiometer.  There  are  two  components  to  the  RSO  IR  signature:  self-emitted  radiation  and 
reflected  IR  radiation.  The  sources  contributing  the  self-emitted  component  of  the  signature  were 
described  in  Equation  (3).  The  reflected  component  of  the  signature  is  predominantly  due  to  reflection  of 
upwelling  Earthshine  from  the  surface  of  the  RSO.  This  is  demonstrated  in  Figure  2.  The  calculation  in 
figure  2  assumes  that  the  Earthshine  is  well  represented  by  a  280  °K  blackbody  for  bands  A  and  C,  and  by 
a  240  °K  blackbody  for  band  D8.  The  reflected  Earthshine  component  also  depends  upon  the  reflectivity 
of  the  RSO  surface  (p  -  1-a)  and  all  of  the  factors  affecting  the  amplitude  of  the  upwelling  Earthshine. 
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Figure  2.  Comparison  of  energy  flux  reflected  from  a  spherical  graybody  RSO  as  a  function  of  altitude  for  three 
sources:  Earthshine  (thin  solid  lines)  and  solar  radiation  reflected  off  the  Earth  ( dashed  lines),  and  direct  solar 
illumination  (thick  solid  lines).  Each  band  is  identified  by  a  letter  and  corresponds  to  the  wavebands  listed  in  Table 
1.  Each  waveband  is  assumed  to  be  square.  RSO  properties  are  listed  in  Table  2.  The  reflectivity  is  l-cq  -  0.5. 


The  relative  importance  of  self-emitted  radiation  and  reflected  Earthshine  can  be  estimated  by 
constructing  a  simple  signature  model.  Under  the  previous  graybody  assumption,  Equations  (2)  and  (3) 
are  solved  for  an  equilibrium  temperature: 


T  = 

eq 


A p  sun  CtFQ  sun  +  FO  Eolith  Kg; 

e*A.a 


—i  1  /  4 


(4) 


which  is  then  used  to  calculate  a  self-emitted  energy  flux  =  oTeq4.  The  upwelling  Earthshine  energy  flux 
at  the  RSO  is  calculated  from  .FOEa[th,  where  the  same  assumption  about  the  temperature  of  the  Earth  used 
to  produce  Figure  2  was  used  to  calculate  Op^h  in  the  three  wavelength  bands.  The  view  factor  can  be 
approximated  as: 


!(. 

An  2 


•cos 


©);©  = 


arcsin 


r  +  h 


(5) 


where  r  is  the  radius  of  the  Earth  (6378  km)  and  h  is  the  RSO  altitude9.  Figure  3  shows  the  result  of  a 
calculation  assuming  a  graybody  spherical  target  and  the  parameters  shown  in  Table  2.  The  plot  shows 
the  ratio  of  self-emitted  radiation  (W/m2)  to  incident  upwelling  Earthshine  (W/m2),  as  a  function  of 
altitude  for  the  SPIRIT  III  A,  C,  and  D  wavebands  (cf.  Table  1).  The  Earthshine  reflected  from  the  RSO 
will  be  a  fraction,  p,  of  the  Earthshine  incident  on  the  RSO;  therefore  Figure  3  represents  an  upper  limit 
on  the  reflected  Earthshine  portion  of  the  RSO  signature.  The  reflected  component  of  the  signature  is 
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inversely  proportional  to  the  distance  between  the  Earth  and  the  RSO.  Therefore,  the  contribution  of  the 
reflected  Earthshine  is  largest  for  RSO  that  are  near  the  Earth. 

The  internal  source  term,  Qmh  was  neglected  in  the  previous  calculation.  This  assumption  is  good 
for  the  reference  spheres,  as  they  have  no  internal  heat  source.  This  assumption  is  probably  adequate  for 
RSOs  in  general.  Debris  and  rocket  bodies  have  no  internal  heat  sources.  Payloads  are  usually  covered 
with  thermal  insulation  blankets  that  facilitate  internal  temperature  control  to  within  ±5  °K,  with  a 
common  internal  temperature  in  the  range  of  288-293  °K10.  The  thermal  insulation  provides  excellent 
thermal  isolation  of  the  interior  of  the  spacecraft  from  the  exterior  environment.  The  amount  of  heat  that 
passes  through  an  insulation  layer  can  be  estimated  by: 

q  -  ZejjfOfT* hot-  T4 coij)  (6) 

where  e#  is  the  effective  emissivity  of  modem  multi-layer  insulation  (about  0.05 10),  Thot  is  the 
temperature  of  the  hot  boundary  of  the  insulation  blanket,  and  Tcoid  is  the  temperature  of  the  cold 
boundary  of  the  blanket11.  If  we  consider  an  aluminized  Kapton  insulation  blanket  with  a  flat  plate  shape 
( (xjeir  -  0.67;  Teq  -  302  °K)12  and  assume  that  the  cold  layer  of  the  blanket  is  at  the  same  temperature  as 
the  interior  of  the  RSO  (-288  °K),  then  the  insulation  actually  lets  external  heat  in.  The  internal  source 
term  does  not  contribute  significantly  to  the  signature  of  the  object.  This  result  is  dependent  upon  the 
exact  optical  properties  of  the  insulation  blankets,  but  provides  a  reasonable  estimate  of  the  transfer  of 
energy  between  the  RSO  interior  and  the  exterior  environment. 

TABLE  2 

Signature  Model  Parameters 
For  a  Spherical  Graybody 


Os 

0.5 

£ir 

6.5 

<*>sun 

1360  W/m2- 

Ao/As 

%  (sphere) 

8 


Altitude  (km) 


Figure  3  Ratio  of  self-emitted  flux  to  incident  up  welling  Earthshine,  as  a  function  of  altitude  and  wavelength  hand 
(A,  C,  and  D)  Each  spectral  band  is  marked  by  its  letter  designation;  A  and  C  overlap  The  RSO  is  assumed  to  be 

a  spherical  graybody . 
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4.  CALCULATION  OF  COLOR  TEMPERATURE  FROM  SPIRIT  III  IR  DATA 


A  common  method  for  calculating  the  temperature  of  a  remote  object  is  to  construct  a  ratio  of  the  flux 
output  from  two  bands  of  a  multi-color  sensor,  and  this  approach  is  taken  with  the  SPIRIT  III  data.  Given 
a  few  assumptions  about  the  nature  of  the  observed  RSO,  this  ratio  can  be  uniquely  related  to  the 
temperature  of  the  target.  The  temperature  so  calculated  is  referred  to  as  the  color  temperature.  The 
procedure  used  is  shown  in  Figure  4.  The  SPIRIT  III  IR  data  arrived  at  the  Surveillance  Data  Analysis 
Center  in  raw  form  (uncorrected  counts).  The  data  were  processed  with  the  SPIRIT  III  CONVERT 
software  in  order  to  transform  the  raw  telemetry  to  calibrated  engineering  units  and  to  identify  and  extract 
the  irradiances  (W/cm2)  and  positions  of  the  observed  point  sources.  A  set  of  post-CONVERT  processing 
filters  identified  which  of  the  point  sources  were  RSOs  and  correlated  them  with  the  known  RSO  catalog. 
This  data  processing  chain  is  known  as  the  Lincoln  Laboratory  SPIRIT  III  Data  Reduction  Pipeline.  The 
pipeline  is  discussed  in  Reference  1.  The  observations  returned  from  the  pipeline  contain  the  RSO  object 
number,  its  position  in  inertial  space  (right  ascension  and  declination),  and  the  observed  in-band 
irradiance  (W/cm2  at  the  instrument)  for  bands  A,  C,  and  D  for  those  point  sources  that  were  correlated 
with  the  RSO  catalog1.  Three  in-band  irradiances  allows  two  band  ratios  to  be  constructed  (A/C  and 
A/D),  in  principle  (cf.,  Fig  4).  However,  the  RSO  signal  in  band  C  was  sometimes  very  weak,  in  which 
case  only  the  A/D-band  ratio  could  be  calculated.  Note  that  we  have  no  a  priori  information  regarding 
the  target’s  temperature;  the  two-dimensional  nature  of  Fig.  4  requires  that  the  points  be  placed 
somewhere.  Second,  the  observed  band  ratio(s)  is  compared  to  a  theoretical  estimate  of  the  band  ratio(s) 
as  a  function  of  the  target’s  temperature.  This  comparison  allows  the  estimated  temperature  of  the  target 
RSO  to  be  quickly  determined  by  interpolation  (cf.  Fig.  4).  We  did  not  design  this  algorithm  to  make  use 
of  the  band  E  data  because  this  waveband  had  a  relatively  poor  calibration1^  and  we  believe  the  resulting 
color  temperatures  would  be  of  poor  quality.  A  future  report  will  address  the  utility  of  calculating  color 
temperatures  using  the  band  E  data. 


Temperature  (°K) 


Figure  4 ,  Observed  band  ratios  A/C  (circle)  and  A/D  (triangle)  compared  with  theoretical  band  ratios  (solid  lines; 
see  text).  Interpolation  of  RSO  temperature  from  the  theoretical  curves  is  represented  by  the  two  arrows.  In  this 
example,  the  A/C  temperature  is  300  °K  and  the  A/D  temperature  is  308  °K. 
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The  assumptions  involved  in  estimating  the  band  ratio(s)  for  the  RSO  target  warrant  discussion. 
The  observed  in-band  irradiance  is  the  spectral  irradiance  of  the  object  integrated  over  the  waveband  of 
interest: 


/,  =4  "kwwmrw 
^  0 


(7) 


where  is  the  irradiance  in  band  /,  R  is  the  range  to  the  RSO,  RSRi(Z)  is  the  relative  spectral  responsivity 
of  band  i  (zero  for  wavelengths  outside  the  band  of  interest),  and  B(XiT)  is  the  spectral  irradiance  of  the 
RSO.  The  relative  spectral  responsivity  incorporates  all  of  the  effects  on  the  RSO  signature  due  to  the 
optics  of  the  observing  instrument  (e.g.,  spectral  responsivity,  transmission  coefficient  of  the  optics,  etc.). 
The  RSR  is  known  for  each  band  due  to  ground  and  on-orbit  calibration14,  but  the  other  quantities  are 
generally  unknown  (R  is  known  only  for  correlated  objects).  Therefore,  it  is  usually  quite  difficult  to 
invert  Equation  (7)  and  obtain  information  about  the  RSO. 

Constructing  a  band  ratio  eliminates  the  dependence  on  Ap  and  R ,  since  those  are  the  same  for 
each  band,  but  leaves  the  dependences  on  the  unknown  emissivity  and  spectral  irradiance  of  the  target. 
Since  the  target  is  assumed  to  be  a  graybody,  its  spectral  irradiance  is  Planckian: 
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where  h  is  Planck’s  constant  (6.63  x  10  M  J-s),  c  is  the  speed  of  light  (2.998  x  10H  m/s),  and  k  is 
Boltzmann’s  constant  (1.38  x  10’23  J/K).  The  graybody  assumption  also  means  that  the  emissivity  is  a 
constant  across  wavebands;  therefore,  when  the  band  ratio  is  formed,  the  emissivity  cancels  out.  The 
remaining  function  implicity  depends  only  on  temperature  and  is  easily  tabulated: 


~\RSRAA)B{A,T)dA 

-  (9) 

\RSRj(A)B(A,T)dA 

0 

The  ratios  RAC  and  Rad  were  calculated  using  these  assumptions,  in  5-degree  increments,  for  the  range  of 
100-3000  °K.  These  ratios  are  compared  to  the  observed  band  ratios  to  derive  the  RSO  color  temperature 
(cf.,  Fig.  4).  Therefore,  the  color  temperatures  are  those  that  the  target  would  have  assuming  that  it  is  a 
graybody,  and  that  all  of  the  observed  signature  is  due  to  the  thermal  self-emission  of  the  target.  The 
consequences  of  these  assumptions  will  be  discussed  later. 

Another  method  of  estimating  the  RSO  temperature  is  to  calculate  the  temperature  of  an  ideal 
blackbody  whose  radiant  energy  is  equivalent  to  the  observed  radiant  energy.  This  is  called  the 
brightness  temperature.  The  brightness  temperature  is  typically  calculated  at  a  specific  wavelength; 
however,  for  our  purposes,  we  define  an  in-band  brightness  temperature  as  the  temperature  of  an  ideal 
blackbody  whose  in-band  radiant  energy  is  equivalent  to  the  observed  in-band  radiant  energy.  If  an  ideal 
blackbody  is  observed  by  a  perfect  sensor  with  square  wavebands,  the  brightness  and  color  temperatures 
will  be  equivalent  and  will  be  equal  to  the  true  temperature  of  the  object.  If  the  sensor  has  non-square 
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wavebands,  then  given  that  we  know  the  spectral  response  of  the  sensor  wavelength  bands  and  that  we 
have  measurements  of  the  radiant  energy  of  a  target  in  those  bands,  the  accuracy  to  which  the  brightness 
temperature  can  be  determined  depends  upon  the  overall  accuracy  of  the  measured  radiant  energy.  The 
accuracy  of  this  measurement  depends  in  turn  upon  the  calibration  of  the  instrument  and  the  noise 
properties  of  the  instrument  and  the  background  scene  against  which  the  target  is  observed.  For  a  sensor 
with  small  measurement  errors,  the  deviation  of  the  brightness  temperature  from  the  true  temperature 
should  be  correspondingly  small.  If  the  target  is  not  a  blackbody  (elr<l),  the  estimated  brightness 
temperature  will  underestimate  the  true  temperature  of  the  object.  For  a  non-black  target  the  color 
temperature  will  be  altered  by  the  ratio  of  ej/e2  (if  e^£2)  where  is  the  graybody  emissivity  in  band  i. 
This  could  result  in  an  under/overestimation  of  the  target  temperature. 
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5.  ACCURACY  OF  THE  COLOR  TEMPERATURE  DETERMINATION 


The  best  possible  test  for  the  temperature  determination  algorithm  is  to  calculate  the  temperature 
of  a  target  for  which  truth  values  are  available.  The  MSX  emissive  reference  sphere  experiments  provide 
such  a  truth  reference.  Five  experiments  wfere  carried  out  over  the  10-month  lifetime  of  the  SPIRIT-III 
instrument;  each  involved  the  deployment  of  a  2-cm  diameter  solid  aluminum  sphere  with  an  emissive, 
Martin  Black  coating.  The  spheres  were  used  to  calibrate  the  SPIRIT  radiometer.  The  spheres  and  their 
deployment  scenarios  were  designed  so  that  their  IR  signatures  would  be  well  characterized,  essentially 
behaving  as  near-ideal  blackbodies  ( o ^  ~  eir  -  0.95).  The  experiments  probed  the  entire  dynamic  range  of 
the  SPIRIT  III  radiometer,  as  the  spheres  moved  away  from  the  MSX  spacecraft.  The  details  on  the 
reference  sphere  experiments  are  summarized  in  Table  3. 


TABLE  3 

MSX  Emissive  Reference  Sphere  Experiment  Parameters 


Sphere 

Deploy  Date 

Initial  Temp  (K) 

Deploy  Velocity  (m/s) 

1 

1996:238:13:34 

259.5 

14.03 

2 

1996:258:11:13 

259.2 

14.40 

3 

1996:327:09:49 

258.9 

13.66 

4 

1996:354:13:05 

258.9 

14.60 

5 

1997:051:13:18 

258.3 

14.20 

The  Data  Certification  and  Technology  Transfer  (DCATT)  PI  group  carried  out  the  reference 
sphere  experiments.  We  obtained  the  observed  in-band  irradiances  from  bands  A,  C,  D,  and  E  and  the 
estimated  MSX-to-sphere  ranges  for  all  5  experiments  from  the  DCATT  PI  group.  The  irradiances  were 
output  of  the  SPIRIT  III  CONVERT  data  reduction  software  package  (Version  6.2.1);  see  Reference  1 
and  the  references  therein  for  more  details  on  the  CONVERT  software  and  processing.  In  addition,  we 
obtained  truth  values  for  the  sphere  temperatures  and  their  in-band  irradiances  from  the  DCATT  group. 
An  example  of  the  reference  sphere  observations  is  shown  in  Figure  5. 

The  truth  values  for  the  spheres’  temperatures  and  irradiances  were  the  result  of  a  significant 
modeling  effort  carried  out  by  MIT  Lincoln  Laboratory1' 1(5  and  the  DCATT  group n  1  \  The  models 
characterized  the  spheres’  thermal  environments,  and  consisted  of  three  components:  incident  solar  flux, 
incident  solar  flux  reflected  from  the  Earth,  and  incident  upwelling  Earthshine  (cf.,  Equation  1).  The 
details  of  the  model  are  covered  in  the  aforementioned  references.  We  note  that  the  DCATT  model  took 
the  variation  of  emissivity  with  wavelength  into  account.  The  1 -sigma  error  in  the  reference  temperature 
model  is  estimated  to  be  <  1.5  K1  . 
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Figure  5.  In-band  irradiance  values  observed  by  SPIRIT  III  during  the  first  emissive  reference  sphere  calibration 

experiment. 


We  used  the  DCATT-supplied  irradiance  observations  and  our  temperature  determination 
algorithm  to  derive  color  temperature  histories,  using  the  A/C  and  A/D  band  ratios,  for  each  of  the 
reference  spheres.  In  each  scene  the  A,  C,  and  D  irradiances  generally  had  slightly  different  time  tags,  so 
the  observed  irradiances  were  interpolated  to  facilitate  formation  of  the  band  ratios  with  a  common  time 
tag.  The  difference  between  the  color  temperature  and  the  truth  model  was  then  calculated  for  each  of  the 
spheres.  The  data  were  screened  for  three-sigma  outliers.  The  mean  difference  and  the  standard 
deviation  of  the  differences  were  calculated;  these  represent  a  bias  and  precision  of  the  temperature 
determination,  respectively. 

Figure  6  shows  the  color  temperature  that  we  calculated  for  the  first  emissive  reference  sphere 
(ERS).  The  triangles  represent  the  A/C  temperature,  and  circles  represent  the  A/D  temperature.  The 
temperature  of  the  sphere  400  seconds  after  deployment  is  about  270  °K,  which  is  near  the  deploy 
temperature  of  259.5  °K.  There  are  no  data  plotted  prior  to  400  seconds,  because  the  pixels  of  the 
radiometer  were  saturated,  rendering  the  irradiance  values  useless.  There  are  three  gaps  in  the  data  due  to 
the  execution  of  routines  designed  to  measure  the  dark  current  of  the  radiometer  focal  planes  (see  also 
Fig.  5).  The  sphere  was  sunlit  when  deployed  but  crossed  the  terminator  at  —  1 600  seconds  after 
deployment;  the  temperature  of  the  sphere  rose  while  sunlit  and  fell  after  the  terminator  crossing,  as 
expected.  The  computed  color  temperatures  are  reasonable;  a  sunlit  sphere  in  thermal  equilibrium,  with 
a=£,  will  have  an  equilibrium  temperature  of  278.5  °K19.  It  is  easy  to  show,  using  the  model  described  in 
Section  3,  that  due  to  the  proximity  of  the  sphere  to  the  Earth,  the  sphere  will  warm  to  an  equilibrium 
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temperature  of  ~287.5  °K;  this  is  similar  to  the  temperature  attained  by  the  sphere  before  it  entered 
eclipse.  The  A/C  and  A/D  color  temperatures  agree  well  with  each  other  during  the  event. 


Time  After  Deployment  (s) 


Figure  6.  " Observed  "  temperatures  for  ERS-I  as  calculated  by  our  algorithm.  The  two  temperature  estimates 
agree  nicely  and  behave  as  expected.  The  error  bars  represent  the  1 -sigma  temperature  error  as  described  in  the 
text  and  in  Table  4.  The  sphere  heats  up  while  illuminated  and  cools  off  after  it  crosses  the  terminator  (-1600  sec 

after  deployment). 


Figure  7  shows  the  three-sigma-screened  temperature  residuals  (TC0|0r  -  T^)  for  ERS-1.  The 
mean  residuals  for  the  A/C  and  A/D  temperatures  are  0.55  and  4.21  °K,  respectively.  There  is  a  small 
positive  bias  in  the  temperature.  The  standard  deviations,  or  precisions,  for  the  A/C  and  A/D  residuals  are 
10.16  and  6.77  °K,  respectively.  The  bias  and  precision  are  added,  in  quadrature,  to  provide  an  estimate 
of  the  total  error  in  the  calculated  temperature  over  the  entire  track: 


Error  =  sj  Bias 2  +  Precision 2  (10) 

For  ERS-1,  the  A/C  and  A/D  temperature  errors  are  10.18  and  7.97  °K,  respectively.  The  biases, 
precisions,  and  errors  for  ERS-1  and  the  other  4  reference  spheres  are  presented  in  Table  4.  A  mean 
accuracy  was  calculated  by  weighting  the  accuracy  for  each  sphere  by  the  number  of  data  points  in  the 
sphere  track.  Averaged  over  all  5  reference  spheres,  the  mean  A/C  (A/D)  temperature  accuracy  is  1 1.02 
(9.32)  °K.  We  believe  that  this  accuracy  conservatively  represents  the  capability  of  the  sensor,  the 
CONVERT  software,  and  the  post-processing  algorithms,  combined.  The  number  of  data  points  used  in 
the  calculation  of  Table  4  and  the  remaining  plots  for  spheres  2-5  appear  in  Appendix  A. 
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Time  After  Deployment  (s) 


Figure  7.  Temperature  residuals  for  ERS-1.  Three-sigma  outliers  have  been  screened  from  the  data. 


TABLE  4 

Emissive  Reference  Sphere  Temperature  Accuracy 


Sphere 

A/C  Bias 
(°K) 

A/C 

Precision 

A/C  Error 
(°K) 

A/D  Bias 
(°K) 

A/D 

Precision 

A/D  Error 
(°K) 

ERS-1 

0.55 

10.16 

10.18 

4.21 

6.77 

7.97 

ERS-2 

3.24 

10.95 

11.42 

7.67 

6.34 

9.95 

ERS-3 

4.17 

10.90 

11.67 

5.07 

10.87 

11.99 

ERS-4 

1.06 

10.23 

10.28 

4.31 

10.09 

10.97 

ERS-5 

-3.77 

11.63 

12.23 

1.50 

8.72 

8.85 
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6.  CALCULATION  AND  ACCURACY  OF  EMISSIVITY-AREA  PRODUCT 


In  addition  to  estimating  the  color  temperature  of  an  RSO,  we  can  use  the  SPIRIT  III  IR  data  to 
estimate  the  emissivity-projected-area  product,  eAp,  for  an  RSO.  An  interesting  technique  for  extracting 
emissivity,  area,  and  temperature  of  a  target,  using  the  IR  signature  observed  by  a  3-band  LWIR  sensor, 
has  been  developed20.  However,  it  is  not  usually  applicable  in  this  work  because  we  do  not  often  have 
observations  in  three  bands.  This  deficiency  is  due  to  the  poor  sensitivity  of  the  SPIRIT  III  C  band  in 
mirror  scan  mode.  Therfore,  we  adopted  a  simple  algorithm,  which  was  employed  to  calculate  eAp  for 
RSOs  detected  by  the  IRAS  satellite12.  The  emissivity-area  may  be  estimated  by  setting  the  in-band 
brightness  temperature  equal  to  the  color  temperature  and  solving  for  eAp21: 


o 


(11) 


In  Equation  (11)  7/  is  the  observed  irradiance  in  band  /,  R  is  the  range  to  the  target,  Tobs  is  the  color 
temperature  calculated  from  the  SPIRIT  observations,  and  B  is  the  Planck  spectral  irradiance  defined  in 
Equation  (8).  We  estimated  eAp  using  the  band  A  data.  This  algorithm  required  that  a  particular  RSO 
observation  be  correlated  to  the  RSO  catalog  (so  that  R  was  known)  and  that  the  RSO  had  been  detected 
in  at  least  two  bands  (so  that  color  temperature  could  be  estimated).  The  ratio  in  Equation  (11)  represents 
the  emissivity-area  product  that  the  target  would  have  assuming  that  it  is  a  graybody  emitter  and  that  all 
of  the  observed  radiation  is  due  to  the  self-emission  of  the  target. 

The  algorithm  for  determining  eAp  was  tested  using  the  ERS  data.  We  assume  the  emissivity- 
area  product  of  the  reference  spheres  is  approximately  equal  to  0.000298  m2.  We  have  assumed  that  the 
area  of  the  2-cm  spheres  is  well  known  and  that  the  emissivity  of  the  spheres  was  0.95r.  Figure  8 
compares  the  eAp  estimates  with  the  truth  value  for  the  first  sphere.  The  estimates  agree  fairly  well  with 
the  truth  value.  We  calculated  bias  and  precision  values  for  each  of  the  spheres  in  a  manner  similar  to  the 
way  we  calculated  the  bias  and  precision  for  the  temperature.  The  biases  and  precisions  were  added  in 
quadrature  to  obtain  an  estimate  of  the  error  in  our  calculation  (cf.,  Equation  9).  The  error  in  the  truth 
model,  mentioned  above,  is  very  small  and  was  neglected.  The  accuracy  estimates  appear  in  Table  5.  For 
clarity,  the  residuals  have  been  normalized  to  the  true  eAp  value  for  the  spheres,  and  have  been  multiplied 
by  100  so  that  the  numbers  in  Table  5  are  percentages.  For  example,  ERS-1  shows  a  bias  of  -9.7, 
showing  that  on  average  our  eAp  estimates  were  ~10%  low  for  this  sphere.  Table  5  shows  that  the 
emissivity-area  estimates  for  the  spheres  are  good  to  only  14-20%.  The  mean  error,  weighted  by  the 
number  of  measurements  per  track,  is  16.8%.  The  eAp  residuals  (calculated  -  truth  values)  for  the  five 
spheres  are  presented  in  Appendix  B. 
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Figure  8.  Estimated  emissivity-area  product  for  ERS-1.  The  solid  black  line  is  the  truth  value  of 0.000298  m  The 
error  bars  represent  the  l -sigma  error  in  the  estimates  as  described  in  the  text  and  Table  5.  These  estimated  were 

made  using  the  SPIRIT  III  A-band  irradiance  data . 


TABLE  5 

Emissive  Reference  Sphere  eAp  Error  Estimates 


Sphere 

eA  Bias 
(%) 

eA 

Precision 

(%) 

eA  Error 
(%) 

ERS-1 

-9.7 

11.6 

15.1 

ERS-2 

-11.2 

9.3 

14.5 

ERS-3 

-6.2 

19.1 

20.1 

ERS-4 

-6.6 

19.2 

20.4 

ERS-5 

+4.5 

16.4 

17.0 

The  precision  of  the  eAp  estimate  and  the  total  error  in  the  estimate  decrease  with  subsequent  sphere 
deployments  (exception:  ERS-5).  Therefore  the  band  ratios  exhibit  increasing  variation  over  the  lifetime 
of  the  radiometer  (this  increased  variation  is  also  seen  by  looking  at  the  color  temperature  errors  as  a 
function  of  time).  This  is  presumably  a  result  of  the  increase  in  the  temperature  of  the  radiometer  focal 
plane  (and  thus  focal  plane  noise)  over  the  lifetime  of  the  instrument.  This  increase  in  focal  plane 
temperature,  and  the  corresponding  increase  in  focal  plane  noise,  was  apparently  accounted  for  by  the 
CONVERT  software  (cf. ,  Ref.  1  and  references  therein)  during  data  reduction.  However,  we  know  of  no 
reasonable  explanation  for  the  increased  variation  in  the  band  ratios  with  time  other  than  an  increase  in 
focal  plane  noise.  It  is  possible  that  the  CONVERT  process  did  not  account  for  all  of  the  variation  in 
observed  irradiance  due  to  focal  plane  temperature  changes. 
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We  also  note  that  the  eAp  estimates  tend  to  be  biased  low  (exception:  ERS-5);  this  is  evident  in  Figure 
8.  If  we  assume  the  reflected  Earthshine  component  of  the  sphere’s  signatures  is  negligible,  which  is  a 
good  assumption  given  the  non-reflective  nature  of  the  spheres,  this  bias  suggets  that  the  estimated  color 
temperatures  were  systematically  high,  or  that  the  MSX-to-sphere  ranges  were  systematically  low  (or 
both).  The  reference  spheres  were  not  tracked  by  the  radars  in  the  Space  Surveillance  Network,  so  we 
have  no  information  regarding  the  actual  time  histories  of  MSX-to-sphere  range  to  compare  with  our 
range  estimates.  Table  4  shows  that  a  positive  bias  in  the  estimated  color  temperatures  does  exist. 
Assuming  that  there  is  no  bias  in  the  temperature  determination  algorithm,  this  temperature  bias  suggests 
that  there  is  a  positive  bias  in  the  band  ratios  or  perhaps  a  uniform  bias  in  the  irradiance  values  that  form 
the  band  ratios.  We  did  not  look  into  the  issue  of  irradiance  errors,  but  the  DCATT  PI  group  has  and  did 
not  observe  a  systematic  bias  in  either  the  band  ratios  or  the  irradiances  for  the  reference  spheres11.  The 
temperature  bias  may  also  be  a  result  of  non-gray  target  behavior.  A  positive  bias  in  the  temperature  can 
indicate  that  £|>£2  rather  than  £i=£2,  where  £1  is  the  spectral  emissivity  in  the  shorter  wavelength  band 
(band  A  in  this  case).  In  the  case  of  the  reference  spheres,  the  graybody  assumption  is  a  good  one  in  the 
6-16  pm  region  of  the  spectrum,  and  the  truth  model  accounted  for  the  non-gray  nature  of  the  Martin 
Black  coating  in  the  other  regions  of  the  spectrum'7  IR.  Therefore  it  seems  most  likely  that  the  bias  in 
temperature,  which  contributes  to  the  bias  in  eA,  is  due  to  an  unresolved  bias  in  the  color  temperature 
calculation  proceedure.  We  note  that  the  weighted  mean  bias  in  the  eA  estimates  is  -5.45%.  This  may  be 
explained  by  an  MSX-to-sphere  range  which  is  systematically  low  by  ~3%,  or  may  be  explained  by  a 
temperature  which  is  systematically  high  by  ~1%.  For  the  ~280  °K  spheres  this  is  about  3  °K  and  is  in 
fair  agreement  with  the  weighted  mean  A/D  color  temperature  bias  of  4.05  °K.  Appendix  C  describes  the 
derivation  of  these  estimates. 
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7.  DISCUSSION 


Two  significant  assumptions  were  made  in  order  to  simplify  the  calculation  of  temperature  and 
eAp  from  the  SPIRIT  III  irradiance  data.  First  was  the  assumption  that  the  target  is  a  graybody.  That 
assumption  is  an  excellent  one  for  the  MSX  Emissive  Reference  Spheres.  In  general,  however,  most 
space  objects  are  probably  not  graybodies.  The  non-gray  characteristics  of  a  target  at  a  specific 
temperature  will  cause  the  irradiances,  and  hence  the  irradiance  ratios,  to  deviate  from  those  of  a  gray 
target  with  the  same  temperature.  The  characterization  of  this  behavior  was  examined  briefly  and  will  be 
expanded  upon  in  the  analysis  of  the  remaining  IR  RSO  data.  These  data  will  be  discussed  in  a  future 
report. 


The  second  assumption  is  that  all  of  the  energy  flux  observed  by  the  SPIRIT  III  radiometer  was 
due  to  the  thermal  self-emission  of  the  target.  This  assumption  simplified  the  calculation  of  theoretical 
band  ratios  (cf.  Section  4)  by  allowing  us  to  neglect  the  contribution  of  reflected  IR  radiation  to  the 
targets’  signature.  The  simple  models  presented  in  Section  3  demonstrate  that  this  assumption  is  a  good 
one  for  targets  with  altitudes  greater  than  ~8000  km  (cf.,  Figure  3),  because  at  these  altitudes,  the  ratio  of 
self-emitted  flux  to  upwelling  Earthshine  is  greater  than  10.  At  lower  altitudes,  the  background  source 
may  contribute  to  the  target  signature  given  that  the  target  is  sufficiently  reflective  in  the  wavelength 
bands  of  interest.  The  reflected  component  will  alter  the  observed  irradiance  ratio  in  such  a  way  that  the 
apparent  temperature  of  the  target  falls  between  the  true  temperature  of  the  target  and  the  temperature 
characteristic  of  the  background  source22.  Techniques  have  been  developed  to  estimate  and  isolate  the 
background  contribution  to  the  signature  and  thus  produce  more  accurate  target  temperature  estimates20. 
Those  techniques  were  not  applied  to  the  ERS  analysis  for  two  reasons:  1 .)  The  spheres  are  not  reflective 
and  the  reflected  component  of  their  signature  should  be  very  small,  and  2.)  The  technique  mentioned 
requires  good  quality  observations  in  at  least  3  different  wavelength  bands.  While  we  have  good 
observations  of  the  spheres  in  bands  A,  C,  D,  and  E,  this  is  not  true  in  general  for  the  other  RSO  targets 
observed  during  the  Surveillance  PI  experiments  which  usually  have  observations  only  in  bands  A  and  D. 
Therefore,  effort  was  not  expended  on  the  implementation  of  these  three-band  algorithms.  Recall  that  we 
did  not  process  band  E  due  to  poor  radiometric  calibration.  We  are  investigating  the  utility  of 
determining  temperatures  using  band  E  data;  should  this  band  prove  useful,  we  will  readdress  the  utility 
of  the  three-band  algorithms. 
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8.  CONCLUSIONS 


This  report  has  described  the  process  used  by  the  SPIRIT  III  data  reduction  pipeline  to  calculate  color 
temperature  and  emissivity-area  product  for  correlated  RSO  observations  obtained  by  the  SPIRIT  III 
radiometer  during  the  surveillance  PI  DCEs.  These  algorithms  were  tested  using  the  best  truth  reference 
available:  the  5  emissive  reference  spheres  deployed  from  the  MSX  over  the  lifetime  of  the  SPIRIT 
instrument.  The  ERS  data  were  fed  into  the  temperature  and  eA  algorithms  to  generate  estimates  of 
those  quantities.  The  estimates  were  compared  to  truth  models  supplied  by  the  DCATT  PI  team.  We 
have  determined  empirically  that  the  1-sigma  temperature  error  for  the  reference  spheres  is  ~1 1°K  for  the 
A/C  band  ratio  and  ~9.3°K  for  the  A/D  band  ratio.  In  addition,  we  have  determined  empirically  that  the 
pipeline  fractional  eA  error  for  the  reference  spheres  is  ~17%.  We  believe  that  these  errors  conservatively 
represent  the  capability  of  the  sensor,  CONVERT  software,  and  post-CONVERT  processing  algorithms. 
In  addition,  the  utility  of  using  band  E  data  is  being  investigated  and  will  be  addressed  in  a  future  report. 
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APPENDIX  A 

EMISSIVE  REFERENCE  SPHERE  TEMPERATURE  PLOTS 


The  color  temperature  calculated  by  the  SPIRIT  III  Data  Reduction  pipeline  for  reference  sphere 
1  (ERS-1)  was  presented  in  Section  5  of  this  report.  This  appendix  presents  the  remaining  reference 
sphere  temperature  data  in  the  same  format  as  Figures  7  and  8.  First  of  all.  Table  A1  presents  the  number 
of  data  points  available  for  the  calculation  of  temperatures  and  emissivity-area  product  (eA). 

TABLE  A1 

Number  of  Data  Points  Used  to  Derive 
Temperature  and  eA  for  the  Reference  Spheres 


Sphere 

Band  Ratio  A/C 

Band  Ratio  A/D 

ERS-1 

201 

220 

ERS-2 

67 

67 

ERS-3 

44 

75 

ERS-4 

40 

73 

ERS-5 

120 

147 

The  top  graph  on  each  of  the  following  pages  is  a  plot  of  the  observed  temperature  as  a  function 
of  time  with  the  A/C  temperatures  represented  by  triangles  and  the  A/D  temperatures  represented  by 
circles.  The  bottom  plot  shows  the  temperature  residuals:  the  observed  temperature  minus  the  DCATT 
truth  temperature.  Three-sigma  outliers  have  been  removed  from  each  data  set.  For  convenient  reference. 
Table  3  is  reproduced  in  Table  A2. 


TABLE  A2 

MSX  Emissive  Reference  Sphere  Experiment  Parameters 


Sphere 

Deploy  Date 

Initial  Temp  (K) 

Deploy  Velocity  (m/s) 

1 

1996:238:13:34 

259.5 

14.03 

2 

1996:258:11:13 

259.2 

14.40 

3 

1996:327:09:49 

258.9 

13.66 

4 

1996:354:13:05 

258.9 

14.60 

5 

1997:051:13:18 

258.3 

14.20 
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Figure  A 1.  Temperature  history  (top)  and  temperature  residuals  (bottom)  for  Emissive  Reference  Sphere  #2.  Bands 
A/C  are  represented  by  triangles  and  circles  represent  bands  A/D. 
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Figure  A2.  Temperature  history  (top)  and  temperature  residuals  (bottom)  for  Emissive  Reference  Sphere  #3.  Bands 
A/C  are  represented  by  triangles  and  circles  represent  bands  A/D. 
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Figure  A3.  Temperature  history  (top)  and  temperature  residuals  (bottom)  for  Emissive  Reference  Sphere  # 4 .  Bands 
A/C  are  represented  by  triangles  and  circles  represent  bands  A/D. 
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Figure  A4.  Temperature  history  (top)  and  temperature  residuals  (bottom)  for  Emissive  Reference  Sphere  #5.  Bands 
A/C  are  represented  by  triangles  and  circles  represent  bands  A/D. 
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APPENDIX  B 

EMISSIVE  REFERENCE  SPHERE  EMISSIVITY-AREA  RESIDUALS 


The  emissivity-area  product  calculated  from  the  band  A  irradiance  by  the  SPIRIT  III  Data 
Reduction  pipeline  for  reference  sphere  1  (ERS-1)  was  presented  in  Section  6  of  this  report.  This 
appendix  presents  the  band-A  eA  residuals  for  all  5  reference  spheres.  Each  plot  shows  the  residual 
(calculated  eA  -  truth  eA)  as  a  function  of  time  after  sphere  deployment.  Three-sigma  outliers  have  been 
removed  from  each  data  set. 


Figure  Bl.  Emissivity-Area  product  residuals  for  ERS-1.  These  data  are  the  differences  between  the  calculated  eA 

and  the  true  eA  for  the  sphere. 
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Figure  B2.  Emissivity-Area  product  residuals  for  ERS-2.  These  data  are  the  differences  between  the  calculated  €A 

and  the  true  eA  for  the  sphere . 


Figure  B3.  Emissivity-Area  product  residuals  for  ERS-3.  These  data  are  the  differences  between  the  calculated  eA 

and  the  true  eA  for  the  sphere. 
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Figure  B4.  Emissivity-Area  product  residuals  for  ERS-4.  These  data  are  the  differences  between  the  calculated  eA 

and  the  true  eA  for  the  sphere. 


Figure  B5.  Emissivity-Area  product  residuals  for  ERS-5.  These  data  are  the  differences  between  the  calculated  eA 

and  the  true  eA  for  the  sphere. 
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APPENDIX  C 

ANALYTIC  ESTIMATE  OF  eA  BIAS 


An  analytic  estimate  of  the  bias  in  eA  due  to  biases  in  range  and  estimated  color  temperature  is 
developed.  We  note  that  the  eA  is  derived  from  the  definition  of  brightness  temperature  as  described  in 
Section  6.  We  repeat  Equation  (11)  here  for  convenience: 


o 


From  simple  propagation  of  errors23,  the  fractional  error  in  eA  can  be  written  as: 
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where  Tc  was  used  for  Tcoior  to  simplify  the  notation,  i  stands  for  wavelength  band,  and  we  assumed  that 
there  is  no  variation  in  the  irradiance  (we  are  interested  only  in  estimating  the  contribution  of  range  and 
temperature  biases  on  eA).  For  range  biases  only,  we  neglect  the  temperature  term  and  note  that  the 
fractional  error  in  eA  is  twice  that  in  range.  For  temperature  bias  only,  we  neglect  the  range  term.  The 
resulting  eA  bias  is  shown  in  Figure  Cl  as  a  function  of  target  temperature  for  a  temperature  bias  of  ATC  = 
+2%.  This  bias  corresponds  to  a  bias  of  -6  °K  at  a  target  temperature  of  280  °K. 


Estimated  Temperature  (=  Tc) 


Figure  Cl.  Estimated  fractional  bias  in  eA  due  to  positive  bias  in  estimated  color  temperature. 
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